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Abstract. The Green function (GF) equation of motion technique for solving 
the effective two-band Hubbard model of high-T c superconductivity in cuprates 
[N.M. Plakida et al., Phys. Rev. B, 51, 16599 (1995); JETP, 97, 331 (2003)] rests 
on the Hubbard operator (HO) algebra. We show that, if we take into account 
the invariance to translations and spin reversal, the HO algebra results in invariance 
properties of several specific correlation functions. The use of these properties allows 
rigorous derivation and simplification of the expressions of the frequency matrix (FM) 
and of the generalized mean field approximation (GMFA) Green functions (GFs) of 
the model. 

For the normal singlet hopping and anomalous exchange pairing correlation 
functions which enter the FM and GMFA-GFs, an approximation procedure based 
on the identification and elimination of exponentially small quantities is described. It 
secures the reduction of the correlation order to GMFA-GF expressions. 
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1. Introduction 

A consistent theoretical model of the high critical temperature superconductivity in 
cuprates is to be able to accommodate both the normal and superconducting states 
under incorporation of the essential features of these systems (see, e.g., p] for a review): 
strong antiferromagnetic (AFM) superexchange interaction inside the Cu02 planes, 
occurrence of two relatively isolated energy bands around the Fermi level, able to develop 
d x 2_ y 2 pairing: one stemming from single particle copper d x 2_ y 2 states and the second 
one from singlet doubly occupied states generated [2] by crystal field interaction; hopping 
conduction for an extremely low density of the free charge carriers. 

The p-d model [3], while incorporating all these features, is too cumbersome and 
cell-cluster perturbation theory [H [5] providing a hierarchy of the various interaction 
terms was used to derive simpler models from it. Extreme limit cases of this reduction 
procedure are various effective one-band t-J models (see, e.g., [SI Ej and references 
therein) which, while unveiling the role played by the AFM exchange interaction in the 
occurrence of the <i-wave pairing, address exclusively the superconducting state. 

The reduction of the p-d model to an effective two-band Hubbard model considered 
by Plakida et al. [8], corroborated with the use of the equation of motion technique 
for thermodynamic Green functions (GF) [9], provided the simplest approach to the 
description of both the normal [H [TO] and the superconducting states [HI [121 02] within 
a frame securing rigorous fulfilment of the Pauli exclusion principle for fermionic states. 

The Green function technique rests on the Hubbard operator algebra. Its rigorous 
implementation onto a system characterized by specific symmetry properties (translation 
invariant two-dimensional spin lattice, spin reversal invariance of the observables) results 
either in characteristic invariance properties of several correlation functions, or in the 
occurrence of some exactly vanishing correlation functions. The use of these results 
allows rigorous derivation and simplification of the expressions of the frequency matrix 
and of the generalized mean field approximation (GMFA) Green functions of the model. 

The obtained expressions contain higher order boson-boson correlation functions 
(CFs). For the CFs involving singlets (normal singlet hopping CFs and anomalous 
exchange pairing CFs), an approximation procedure which avoids the usual decoupling 
schemes and, yet, secures the correlation order reduction to GMFA-GF expressions, 
under the identification and elimination of exponentially small quantities, is described. 

The organization of the paper is as follows. Sec. [2] summarizes essentials of the 
two-band Hubbard model and GMFA-GF equations. Sec. [3] describes the invariance 
properties following from the translation invariance of the underlying spin lattice. 
Sec. S] derives invariance properties and constraints following from the invariance of the 
macroscopic properties of the system under spin reversal. On the basis of the results of 
Sec. [3] and HI rigorous derivation of the frequency matrix in the (r, ^-representation is 
done in Sec. [51 The derivation of GMFA-GF expressions for the boson-boson correlation 
functions involving singlets is discussed in Sec. El 

Collecting together the results of sections [5] and [SI expressions of the frequency 
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matrix and of the GMFA Green function matrix are derived in the (q, (^-representation 
in sections [7] and [8] respectively. These results explicitly incorporate both hole-doping 
and electron-doping features of the cuprate systems through the singlet hopping and 
superconducting pairing terms. 

The paper ends with conclusions in section [9j 

2. Mean field approximation 

The Hamiltonian of the effective two-band singlet-hole Hubbard model [8] is written in 
the form 

+ )c n 2^T 1 / +/c 22 2^ r i/ + ^21 2^ 2(T ( r i/ + t m ) C 1 ) 

i,a i,a i,<r 

The summation label i runs over the sites of an infinite two-dimensional (2D) square 
array the lattice constants of which, defined by the underlying single crystal 

structure. The spin projection values in the sums over a are a = ±l/2,a = —a. 

The Hubbard operators (HOs) X^ = \ia)(i/3\ are defined for the four states of the 
model at each lattice site i: |0) (vacuum), \a) = |f) and \a) = \i) (single particle spin 
states inside the hole subband), and |2) = |||) (singlet state in the singlet subband). 

The multiplication rule holds Xf X] v = Sp^Xf 71 . The HOs may be fermionic (single 
spin state creation/annihilation in a subband) or bosonic (singlet creation/annihilation, 
spin or charge densities, particle numbers). For a pair of fermionic HOs, the 
anticommutator rule holds {X™ , Xj V } = 5ij{5^Xf v + 8 va X^) whereas, if one or both 
HOs are bosonic, the commutation rule holds [Xf 13 , Xj v ] = dij^^Xf 11 — 5 m X^). At 
each lattice site i, the constraint of no double occupancy of any quantum state \ia) is 
rigorously fulfilled due to the completeness relation X®° + X° a + X° a + Xf 2 = 1. 

In (Op) , Ex — en — \i denotes the hole subband energy for the renormalized energy id 
of a (i-hole and the chemical potential /i. The energy parameter of the singlet subband 
is E 2 = 2E 1 + A, where A w A pd = e p — Ed is an effective Coulomb energy U e ff 
corresponding to the difference between the two energy levels of the model. 

In the description of the hopping processes, the label 1 points to the hole subband 
and 2 to the singlet subband. The hopping energy parameter K a b = 2t p dK a b depends on 
t pd , the hopping p-d integral, and on energy band dependent form factors, K ab . Inband 
(/Cn,/C 2 2) and interband (/C 2 i = /Ci 2 ) processess are present. The Hubbard 1-forms 

r ^ = Y^ VimX fX^ (2) 

incorporate the overall effects of specific hopping processes (through the labels (a/?, 777) 
of the pair of Hubbard operators) involving the lattice site i and its neighbouring sites. 

Up to three coordination spheres around the reference site i do contribute [U [5] to 
the sum (T5]), each being characterized by a small specific value of the overlap coefficients 
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Uij {yi for the nearest neighbour (nn), v 2 for the next nearest neighbour (nnn), u 3 for 
the third coordination spheres). 

The quasi-particle spectrum and superconducting pairing for the Hamiltonian (DO) 
are obtained [TTJ [12] from the two-time 4x4 GF matrix (in Zubarev notation [9]) 

G ija (t - t') = ((X ia (t)\xl(t'))) = -i9(t - t')({X la (t),Xl}), (3) 

where (■ • •) denotes the statistical average over the Gibbs grand canonical ensemble. 
The GF 03]) is defined for the four-component Nambu column operator 

X m = {X? Xf Xf XfY (4) 

where the superscript T denotes the transposition. In ([3]), = (Xf X?° X" 2 Xf) 
is the adjoint operator of Xj a . 

The GF matrix in (r, to) -represent at ion is related to the expression ([3]) of the GF 
matrix in (r, t) -represent at ion by the non-unitary Fourier transform, 

+oo 

G lJa (t ~t') = ^ J G ijcT (u) e-^W duj . (5) 

— oo 

The energy spectrum of the translation invariant spin lattice of ([T]) is solved in 
the reciprocal space. The GF matrix in this (q, ^-representation is related to the GF 
matrix in (r, (^-representation by the non-unitary discrete Fourier transform 

^H = ^E e " iq(rj " r!) ^( ( i^)- ( 6 ) 



N 

q 



For an elemental GF of labels (a/3, 77/), we use the notation ((X® 13 (t)\Xf (t'))) in 



the (r, t) -represent at ion and, similarly, ((X° (assuming Hubbard operators at 

t = 0), in the (r, (^-representation. In the (q, (^-representation, it is convenient to use 
the notation G a ^{q, u). 

We shall consider henceforth the GMFA-GF, G°(q, uj). Its derivation involves: 

(i) Differentiation of the GF ([3]) with respect to t and use of the equations of motion 
for the Heisenberg operators X^ (t). 

(ii) Derivation of an algebraic equation for Gij a (u>), Eq. ([5]). 

(iii) Elimination of the contribution of the inelastic processes to the commutator 
Z ia = [X ia , H] entering the equation of motion of Gij a {uj). 

(iv) Transformation to (q, (^-representation of the obtained equation of G^ a (uj) by 
means of the Fourier transform (J6]). 

This finally yields 

Gl (q, u) = x [xw ~ X(q)] \ , (7) 
X={{X ia ,Xl}), (8) 

• / ^cr(q) ^ 6 ^ ^ J ^ vAjjCT, T{j Tj Fj , (9) 
r u 

A ll(7 = ({[X la ,H},Xl}). (10) 
The matrix Aij a is Hermitian. 
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3. Translation invariance of the spin lattice 

Four consequences follow from the translation invariance of the spin lattice. 

• The definition of the Hubbard 1-form ([2]) over a translation invariant spin lattice 
results in the identity (which secures the hermiticity of the Hamiltonian H): 



• The Green function ([3]) of the model Hamiltonian (PJ) depends only on the distance 
r ij = \ r j~ r i\ between the position vectors at the lattice sites % and j [9]. 

• The one-site statistical averages are independent on the site label i, (X?) = (X^), 
(V For this reason, the site label in the one-site averages will be omitted. 

• The two-site statistical averages {Xf Xj V ) remain invariant under the interchange 
of the site labels i and j, 

(XfX]V) = (XfXD, i^j (12) 
4. Spin reversal invariance 

The energy spectrum of the system described by the Hamiltonian (CQ) does not depend 
on the specific values a = ±1/2 of the spin projection. As a consequence, the definition 
of the GF ([3]) either in terms of the cr-Nambu operator (J4j) or the cx-Nambu operator 

X- a = (Xf X? Af°) T (13) 

has to result in mathematically equivalent descriptions of the observables. This means, 
however, that the mathematical structures of the frequency matrices Aij a , Eq. (TlOl) . and 
Aija = ({[Xia, H], Xj-}) emerging from the cx-Nambu operator ({TBI , have to be related 
to each other. 

The identification of the existing relationships is constructive: we calculate and 
compare the corresponding matrix elements of Aij a and A^a- The multiplication rules 
and the commutation/anticommutation relations satisfied by the Hubbard operators 
result in the following general expression of the elemental anticommutators entering 
their definitions: 

{[X* H],X^} = 5 l3 C^ + (1 - S^jT^, (14) 



with one-site contributions given by 



X^ + 



a 

+5 X2 (-E 2 Xf* + /C 22 <f° + K> 21 ^r^°) - 

cr a 



+ 
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VtT,OtT\ 



i \ V r 771 \rua- 1 tr i/Oflcr i^- 1/2, 2a , 1^- r> I u2,a0 , v0,o2 

+ 2 J V[ £ i I i + /c n T "i,i -^r^j +/C21 ■ 2o-(r M +r M 

+ /c 22 £ <r 2 + /c 21 x; 

cr cr 

-w^n £ <r° + /c 21 J2 JK? 5 *)} 



+ 



H T l,j 


+ 2a/C 21 rr;- 25 ) 






n r i,i 


+ 2a/C 21 r 1 i ;r 2 ) 


- 'W^T^' 2 ' 7 


- 2^^°) 


A<^,0o- 
H T l,i 


+ 2a/C 21 r^ 2 ) 


c / 1^ Atp,2cr 


- 2aJC 21 T*?> B0 ) 




+ 2a/C 21 r 1 A f 2 ") 




- 2a£ 21 r#'°*) 



+ 



and two-site contributions given by 



+ 



+<W [(^5A ff )(-/c 11 x^x;°+/c 22 x^xf)+(5 A0 /c 11 -5 A2 /c 22 )^xrx 

cr 



va 
3 



+JC 2 v2a 



s^x^xf+s^x^xf+s^x^x^+x^x]^ } + 



+ ^cr ^•pO^llX^X^ — S tp 2K-22X^ 2 Xj' 7 + 



+/C 21 -2a 



cv XfXf + c^X^Xf + 6 v ^X^Xf + X™X. 



• AO t^ ""/' 



A2 i v-AO v^2\ 



+ 



+y]S ucr UxoJCnX^ XY -byaK-nX^Xf +K, 2 \-2a(5^XY Xf +5 Al2 X l ACT X^ 
~y^<W ^ £iiX Act XJ°-<^ 2 /\; 22 X Act ^^ 

cr 

+^2i 2cr {dxodrtX? 11 XJ* 3 - 5x2$uqX^ X^ - S^d^Xf'X™ + d^d^X^X. 



v0\ 



+ 



The comparison of the results obtained from (fill) for the corresponding matrix 
elements of Aij a and Aija and the use of the translation invariance properties (fTTT) 
and ( fl2l) result in four distinct kinds of relationships: 



Under the spin reversal a — > cr, the following invariance properties hold for the 
normal one-site statistical averages: 



<*r> = w> 



(15) 



Mean field Green functions of Hubbard model of superconductivity 



7 



2a(T%*°) = 2a(f f °> (17) 



The identity (C° ' a +C i ) = holds, therefrom we get for the one-site anomalous 
averages, 

(X° 2 > = (18) 

«r 2 > = -«r 2 > (19) 

(r^ *) = «?' ff2 > (20) 

The first two equations imply that the contributions of the one-site terms (Xf 2 ) 
and ^2 tT (Ti a i ' a2 ) to the superconducting pairing vanish identically irrespective of the 
model details (like, e.g., the relationship between the lattice constants a x and a y ). 
For a rectangular spin lattice (a x ^ a y ), Eq. (120]) points to the occurrence of a small 
non-vanishing one-site contribution to the superconducting pairing originating 
equally in both energy subbands. However, over the square spin lattice (pQ) 
(a x = a y ), each term of (1201) vanishes for <i-wave pairing due to the symmetry 
in the reciprocal space [T2|. 

Under the spin reversal a — > a, the following invariance properties hold for the 
two-site statistical averages: 

(X^Xf ) = {X'fX'f ), (Xl a Xf) = {XfX™) (21) 
(XfXY) = (Xf 2 Xf ), (X 4 00 XJ CT ) = {XfXf ) (22) 
(Xf X™) = (Xf Xf ). (23) 

The operator of the number of particles at site i within the singlet subband, N iy is 
the sum of spin a and a components, 

N^n^ + n^, n iCT = Xf + Xf, n i& = Xf + X 22 . (24) 

Similar relationships hold for the number of particles at site i within the hole 
subband, X/ 1 , 

N? = nl + n%, nt = Xr + X°°, n% = Xf + X°°. (25) 

Due to the completeness relation, 

Ni + = 2, n m + nt = n id + = 1. (26) 

These equalities simply reflect the fact that, at a given lattice site i, there is a single 
spin state of predefined spin projection, whereas the total number of spin states 
equals two. 

Therefore, the operator Ni, Eq. (124"1) . provides unique characterization of the 
occupied states within the model [8l [T2| [TO]. 
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5. Frequency matrix in (r, ^-representation 

A straightforward consequence of the results established in section H] is the simplest 
general expression of the frequency matrix Aij a , Eq. ffTUj) : 

« -Ivh'Mifcy -fc T )' (27) 

The one-site 2x2 matrix c a is Hermitian, its elements do not depend on the particular 
lattice site i, 



(2f 



(Ei + A)x 2 + «22 2cra 2 i 

and are expressed in terms of the spin reversal invariant quantities 

X2 = (n ia ) = (riia) (29) 

Xi = «> = (t&> = 1 - xa (30) 

/ 0<r,crO\ i^- / a2.2a\ /oi\ 

o 2 2 = AJii(r 1 )- ^22(71 ) (31) 

Q21 = (/C n - /C 22 ) • 2a(rr 2 ^ ) + /C 21 ((rf • 50 > - (rf • 2 < J ». (32) 
The normal hopping 2x2 matrix Z)j 3£T is symmetric, 
( d 22 2ad 21 " 

D '- = t 1 (33) 



Due to the constraints (1211) -( 1221) . the charge-spin correlations entering the matrix 
elements of f )33|) get exactly decoupled from each other, such that 

C = K 22 (x'+x£)-£iix^ 



d}} = lC n \xij + (Xi - X2H + xj] - ^22Xtj h 



df; = /c 21 [(xii - X2^0 + xf;] - Kaa'i > 

with the three spin reversal invariant weighted boson-boson correlation functions 

representing respectively charge-charge (c), spin-spin (S), and singlet-hopping (s-h) 
correlations: 

X % =v ij {N i N J )/4, (34) 

Xij = Vij{SiSj) (35) 

X? = MX?*?) ( 36 ) 
In (E5J), S, = (5?, Sf), with S* = (Xr - Xp/2 and Sf = Xf. 

The anomalous hopping 2x2 matrix Ay CT has a very special form namely, 

x _ I — /C 2 i ■ 2a |(/Cn + /C 22 ) \ pa j r , . 

\ -|(/C n + /C 22 ) " /C 21 -2a ) Xi i [6() 

where the spin reversal invariant weighted boson-boson pairing (pair) correlation 
function is given by 

XT = = toviXfiXr + X f)) = ( 38 ) 

= - vv(N}X?) = -2u v ({X?> + Xf)X™ ). (39) 
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In Eqs. f)38p and (1391) . the derivation of the second expression from the first one makes 
use of the spin reversal invariance property (123]) . 

To get a workable expression of the frequency matrix, approximations have to be 
derived for the boson-boson statistical averages entering the two-site hopping matrix 
elements. In the next section we show that the method of reference p2], yielding 
the pairing correlation function (X® 2 Nj) in terms of GMFA Green functions within 
an approach able to identify and rule out exponentially small terms, can be extended to 
the singlet hopping correlations (X^X 20 ) as well. 



6. Hopping processes involving singlets 



The right approach to the reduction of the order of correlation of the boson-boson 



statistical averages (Xf 2 X^) = (X^Xf 2 ) goes differently for the hole-doped and 
electron-doped cuprates. 

• Reduction of the correlation order for hole-doped cuprates 

In these systems, the Fermi level (the zero point energy) stays in the singlet 
subband. We get the estimates E 2 ~ -A, E 2 - A ~ -2A, E 2 + A ~ 0. With 
A ~ 3eV, /3A ~ 3.5 • 10 4 T _1 . Therefore, at T < 300K, the quantities containing the 
factor e^ 2 ~ e _/3A < e~ 100 < 1CT 44 are negligible. 

We start with the following form of the spectral theorem [9] 



I y02 _ 1 f da; 



(40) 



written for anticommutator retarded (u + ie), respectively advanced (u — ie) Green 
functions. Their equation of motion in the (r, ^-representation is 

{u-E 2 ){{X™\X^ » w ~ 2{Xf X^)+/C 21 ^2a[((r 1 °f CT |X^)) W - ((r^ 2 |X" M )) J (41) 



where, for the sake of simplicity, the labels ±ie,e = + , describing respectively the 
retarded and the advanced Green functions have been omitted. In Eq. fj4T|) . the higher 
order r.h.s. contributions coming from the inband hopping terms have been dropped off. 
Replacing (T4T]) in (HOI , we get 

da; 



/+oo 
T 
-oo 



x [ — — )lm 



+ e"^ 
1 



Oct, Oct i x^A/j 



-to — E 2 + ie 

To evaluate the imaginary part, we use the identity [9] 

1 _ 1 



uj — E 2 + ie 



V 



uj - En 



— m5{uj — E 2 



The integrals over the ^-function yield (finite) GF real parts at u = E 2 , multiplied by 
a thermodynamic factor ~ e _/3A <C 1. The imaginary part of the hole subband GF 



X£ IP ~ (1 - *<,-)— ■ 2<r(Xf Xf >, (44) 
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(( r il' ° \ X j^)) w+ie shows a 5-like maximum at uo = E 2 — A, where {ui — E<i)~ x ~ A -1 and 
the thermodynamic factor reaches a value ~ e~ 2A . The only non- negligible contribution 
to the principal part integral comes from the singlet subband GF ((r^' <j2 |X J AM )) tl , + i e the 
imaginary part of which shows a 5-like maximum at u = E 2 + A ~ 0. This allows us to 
approximate (u — E^ 1 ~ A -1 within the integral over the singlet subband GF to get 

(x?x») ~ (1 - E ^(nr x h (42) 

Replacing this result in Eq. and using (J2D we get 

X ff ~ (1 - <%) \Au iy 2a(XfXp) - £ 2^f 2 X* 2 X J >] (43) 

Omitting the three-site terms, we get the two-site approximation of the superconducting 
pairing originating in the singlet subband, 

4/C 21 4 

which reproduces the well-known two-site exchange term of the t-J model. 

For the singlet hopping correlation function, (j42p yields the two-site approximation 

Xl~ h ~ (1 - 5,,)^4 • 2a(Xf 2 Xf > (45) 

• Reduction of the correlation order for electron-doped cuprates 

The Fermi level (the zero point energy) stays now in the hole subband. We have 
the estimates E 2 ~ A, E 2 + A ~ 2A, E 2 - A ~ 0. 

It is convenient now to start with the alternative form of the spectral theorem [9] 

+oo 

= ^ / [<<*W> W ~ ((XFIX*))^] , (46) 

— oo 

with the retarded and advanced GFs following from the same equation (j4Tj) . 

Exponentially small quantities result from the 5-term of (u — E 2 + is)" 1 and from 
the singlet subband GF ((r^ 2,<j2 |X J AM )) w+ie . The hole subband GF ((r[ l J' 0f7 |X J A/i )) w+ i e , 
yields the non-negligible contribution 

(X*X?) ~ (1 - <%) ^ E 2^(^ M <f CT ) (47) 

Replacing in (1391) and omitting the three-site terms, we get the two-site approximation 
of the superconducting pairing originating in the hole subband, 

A fr"" 2 

XT * (1 - <%) ■ 2<r(Xf Xf ) (48) 
Finally, the two-site approximation of the singlet-hopping correlation function is 

(Xf Xf ) * (1 - ^)?^ • 2a(Xf Xf ). (49) 

In conclusion, the GMFA superconducting pairing is a second order effect. The 
lowest order contribution to it originates in interband hopping correlating annihilation 
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(or creation) of pairs of spins at neighbouring lattice sites % and j within that energy 
subband which crosses the Fermi level. 

Similarly, the singlet hopping is a second order effect as well. It mainly proceeds 
by interband % T^- j single particle jumps from the upper energy subband to the lower 
energy subband. 



7. Frequency matrix in (q, ^-representation 

The calculation of the matrix elements of Ar(q) from Eq. ([9]) asks for three essentially 
different kinds of Fourier transforms, namely, 

• The averages of the Hubbard 1-forms entering Eqs. ( 13T|) and ( [321) result in sums of 
products of q-space averages and geometrical form factors: 

(rtm = E "« ■ JJ £<* A ^> q 7 Q (q) (50) 

a=l q 

for label sets {(A/i, uip)} e {{0a, a0); (a2, 2a); (cx2, a0)}. 

The quantity (X A ^X^) q denotes the average of the q-space image of the product 
of Hubbard operators of labels A/i and vtp respectively, 

{X x, X u V) = J_ f uj + ie)- G x ^(q, u - is)] (51) 

Z7T J 1 + L J 

— CO 

Finally, in Eq. ([50]) . 7 Q (q) denote the nn (a = 1), nnn (a = 2), and third 
neighbour (a = 3) geometrical form factors, 7i(q) = 2[cos(q x a x ) + cos(q y a y )}, 
72 (q) = 4:COs(q x a x ) cos(q y a y ), 73 (q) = 2[cos(2q x a x ) + cos(2q y a y )]. 

• For the two-site weighted singlet hopping (136|) and the superconducting pairing (138|) . 
the Fourier transforms result in convolutions of specific averages and geometrical 
form factors. The results are as follows: 

— Singlet hopping 

X S " h (q) - E < -jf E S ^7a(q - k) (52) 

a=l k 

where S k = 2a(X a2 X°°) k , while S k = 2a(X°°X 2a ) k for hole-doped and electron- 
doped cuprates respectively, with averages defined in ( 15T1) . 

— Superconducting pairing 

X pair (q) = E ^ -jf E n ^(q " k) (53) 

a=l k 

where n k = 2a(X a2 X° 2 ) k , while n k = 2(x(X 0s X 0<T ) k for hole-doped and electron- 
doped cuprates respectively, with averages defined in (15TT) . 

• The charge-charge and spin-spin correlation functions (|34p and fl35l) are treated 
approximately following [H [10] : 
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- The order of the charge-charge correlation function (JVjiVj-) is lowered using a 
Hubbard type I approximation decoupling procedure (NiNj) ~ (Ni)(Nj) = 2% 2 - 

- The spin-spin correlation function (S»S 3 -) is kept undecoupled, but treated 
phenomenologically. Eq. ([2]) implies the occurrence of up to three non-vanishing 
spin-spin correlation functions: nn, xf = (SiS i±ax/y ), nnn, xf = (SiS i±ax ± ay ), and 
X§ = (SiSi±2a x/y ) . These are site independent quantities. 

Using the above results, we get from (Q and (1271) the mathematical structure of the 
frequency matrix Ar(q) as follows, 

A(a)=( ^ (q) ^ (q) ) (54) 
The normal 2x2 matrix contributions to A a (q) show the characteristic cr-dependence, 

P I \ I C 22 2crC 21 \ /A / xxT I _C 22 2aC* 21 \ , . 

E ^ = { 2a 4l c u J ; = ( 20^ -c„ J (55) 

with the a-independent terms c a fe carrying normal one-site and two-site matrix elements, 

c 22 = c 22 (q) = (#i + A)x 2 + a 22 + d 22 (q) 
en = cn(q) = #iXi + a 22 + dn(q) 
c 2 i = c 2i (q) = a 2i + d 2 i(q) 

3 1 

d ab (q) = JC ab J2^1a(q)bC S a + (-l) a+b XaXb} + ^-W^O) 
a=l 

The one-site terms are defined by Eqs. ( l31~T) - (!32l) and ( J50l) . The exchange energy 
parameters are given by 

J ab = 4/C ab /C 21 /A, {ab} G {22, 11, 21}, (56) 

while the singlet hopping contribution x s ~ h (q) is given by Eq. fl52l) . 

The anomalous 2x2 matrix contributions to ^ CT (q), obtained from (|37|) . show the 
characteristic a-dependence, 

* (q) = ( 66 ^ ■ (* (q)) f = ( ^ (5T) 

^ 2^1^' 1 fflqJJ ^ 66* 2a^b* J [ > 

with £i = J 2 i, £ 2 = (^11 + J22V2, whereas 6 = o(q) is a shorthand notation for the 
pairing matrix element fl53l) . 

Remark 1 T/ie spin reversal a — > er symmetry properties of the elemental Green 
functions entering the matrix GF are identical to i/iose established for the underlying 
frequency matrix A a (q). 
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From Eqs. (1151) and (1181) it follows that the matrix x, Eq. (jSJ), is diagonal and spin 
reversal invariant, with two nonvanishing matrix elements, 

\= ( i °, )■ v = (:;")■ 0= ( ;; ;; )• <™> 





o X 

where X2 and Xi are given by Eqs. f[2"9l and (13"0"j) respectively. 

Replacing in ([7]) the expressions (|58|) of the matrix x an d (|54|) of the frequency 
matrix ^4 CT (q), we get a structure of the GMFA-GF matrix obeying the general symmetry 
properties established in [IT] . 



\ raq,^)) T -(G>(q,-^)) 

where the argument u; carries, in fact, the complex value uo + ie,e = + . (Hence the 
elemental GFs containing the argument uo point to retarded GFs, while those containing 
the argument —uo point to advanced GFs.) 

The normal 2x2 matrix G>(q, to) shows the characteristic a-dependence, 

with the a-independent components g ab (q,uo) found from 

g ab (q, uo) = A ab uo 3 + B ab uo 2 + C ab uo + D ab , {ab} G {22, 11, 21}. 
Here the coefficients A ab are given respectively by 

A22 = X2, A n = xi, Mi = 0, 
while B ab , Cab, D ab are q-dependent coefficients: 

5 22 (q) = c 22 , Bn(q) = cu, S 2 i(q) = c 21 

C 22 (q) = " [XsCc?! + C?l&| 2 ) + Xi(|c 21 | 2 + C||6| 2 )]/X? 
Cu(q) = - [Xi(4 + £ 2 |&| 2 ) + X 2 (|c 21 | 2 + £ 2 2 |&| 2 )]/^ 

C 2 i(q) = [c 2 i(x 2 cn + xic 22 ) - 66l^l 2 ]/(xix 2 ) 

A 2 (q) = -[cii(c 22 c n - |c 21 | 2 ) + (c 22 £ 2 + c n il + 23?(c 21 )66)l&| 2 ]/Xi 
Ai(q) = -[c 22 (c 22 c n - |c 21 | 2 ) + (c n £ 2 + c 22 £ 2 2 + 2^(c 21 )£ 1 £ 2 )|&| 2 ]/ X 2 . 
£>2i(q) = {c 2 i(c 22 c n - |c 21 | 2 ) - [c*^ 2 + c 21 £ 2 2 + (c 22 + c n )ei6]|&| 2 }/(XiX 2 ) 
The anomalous 2x2 matrix F°(q, c<j) shows the characteristic a-dependence, 

F°(n,A=( 2(T /2 2 (q,^) / 2 i(q,w) ^ 1 

J l -/ 2 i(q,-^) 2a/ u (q,w) i'^(q,a;) 
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with the elemental GFs / b(q, oS) given by 

f aa (q,u) = (P aa u 2 + R aa )b, {aa} G {22, 11}, 
/ 21 (q, u) = {P 21 u 2 + Q 2l u + R 21 )b. 

Here, P 22 = — £i, -Pn = £i, and P 2 i = —£2 are q- independent, while 

i? 22 (q) = [(cL + 4)^1 + 2 CllC21 e 2 + 6(^1 - iim'Vxi 

^u(q) = -[(4 + 4i 2 )6 + 2c 22 c 21 6 + £i(£ 2 - el)l&l 2 ]/xl 

i?2i(q) = [(cue*! + c 22 c 21 )£i + (c 22 c n + |c 21 |% - £ 2 (£ 2 - il)\b\ 2 } / {X1X2) 

<?2i(q) = [(X2C21 - Xic 2 i)6 + (X2C11 - Xic 22 )6]/(XiX 2 )- 

The denominator £>(q, u;) occurring in Eqs. ( l60l) and ( 16TI) . which is proportional to the 
determinant of the matrix x^ ~ Ar(q) i n ©> shows the following monic bi-quadratic 
dependence in u: 

X>(q, u;) = (cu 2 — utu + f )(w 2 + mcj + t>), (62) 
where t> = u(q) and u = u(q) are found respectively from 

v 2 = { [(c 22 c n -|c 21 | 2 ) - (£ 2 - e 2 )\b\ 2 ] 2 + [[(c 22 + c n ) + 23*(c 21 )] 2 £ 2 - 

-4(c 22+Cll )3?(c 21 )6(6 - 6) " 4|c 21 | 2 te 2 -e 2 2 )] |&| 2 }/(x?xl) (63) 

^ 2 - 2^ = + e,\b\ 2 ) + ^(c 2 2 + £ 2 |6| 2 ) + — (|c 21 | 2 + £ 2 |6| 2 ). (64) 

xi xi X1X2 

A necessary consistency condition to be satisfied by the parameters of the model at any 
vector q inside the Brillouin zone is v 2 (q) > 0. 

Remark 2 The zeros of the determinant of the GMFA-GF, 

Z>(q,w) = (65) 

provide the GMFA energy spectrum of the system. 

At every wave vector q inside the Brillouin zone, this yields for the superconducting 
state the energy eigenvalue set 

{^(q), fi 2 (q), -n 2 (q), -fli( q)}, 

fi 1)2 (q) = («/2) ± v 7 ^) 2 -^. (66) 

In the normal state (6 = 0), Eqs. fl63|) and ([51]) reduce respectively to 

^0 = (c 22 /x 2 )(cii/xi) - |c 2 i| 2 /(xiX 2 ) 

Mo = (c 22 /x 2 ) + (cn/xi) 

such that the energy spectrum is given by the roots of the second order equation 
oj 2 — Uqoj + v = solved in [8] . 

Finally, if we assume a pure Hubbard model (i.e., energy band independent hopping 
parameters, fCn = /C 22 = /C 2 i = t, [10]), then a significant simplification of the 
equations derived in the last two sections is obtained. The normal 2x2 matrix E a (q) 
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becomes symmetric and so is the normal GMFA-GF G°(q, 10). Moreover, there is a 
single exchange energy parameter in (|5T1) . £1 = £2 = J — 4i 2 /A, which simplifies the 

anomalous 2x2 frequency matrix to $ CT (q) = \ ^° ] J6, such that the quantities 

u and v in the expression (162D of the GF determinant reduce to 

^ 2 = [(c 22 c n - 4) 2 + (c 22 + c n + 2c 21 )V 2 |6| 2 ]/(x?xl) (67) 
M 2 - 2v = ixlcj, + x?4 + 2xiX 2 c| + J 2 |6| 2 ]/(X?^)- (68) 

A non-negative value t> > always follows from Eq. (157)1 . however, the reality of the 
solutions (1661) needs investigation of the domain of variation of the adjustable parameters 
of the model. 



9. Conclusions 

The two-band Hubbard model of the high T c superconductivity in cuprates [H EE] uses 
Hubbard operator algebra on a physical system characterized by specific invariance 
symmetries with respect to translations and spin reversal. 

In the present paper we have shown that the system symmetries result either in 
invariance properties or exact vanishing of several characteristic statistical averages. The 
vanishing of the one-site anomalous matrix elements is shown to be a property which 
is embedded in the Hubbard operator algebra. Another worth mentioning consequence 
following from the spin reversal invariance properties of the two-site statistical averages 
is the exact decoupling from each other of the charge and spin correlations entering 
the matrix elements of the frequency matrix. The use of these results allowed rigorous 
derivation and simplification of the expression of the frequency matrix of the generalized 
mean field approximation (GMFA) Green function (GF) matrix of the model. 

For the higher order boson-boson averages (A? 2 A 20 ) and (Xf 2 Nj), which enter 
respectively the normal singlet hopping and anomalous exchange pairing contributions 
to the frequency matrix, an approximation procedure resulting in GMFA-GF expressions 
was described. The procedure avoids the current decoupling schemes [El fT5] . Its 
principle, first formulated in [T2J, consists in the identification and elimination of 
exponentially small contributions to the spectral theorem representations of these 
statistical averages. 

A point worth noting is that the proper identification of exponentially small 
quantities asks for the use of different starting expressions of the spectral theorem for 
the hole-doped and electron-doped cuprates. 

The results of the reduction procedure may be summarized as follows: 

• The singlet hopping is a second order effect which may be described as interband 
iT^- j single particle jumps from the upper to the lower energy subband. 

• The GMFA superconducting pairing is a second order effect, the lowest order 
contribution to which originates in interband hopping correlating the annihilation 
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(creation) of spin pairs at neighbouring lattice sites % and j within that energy 
subband which crosses the Fermi level. 

The derivation of the most general and simplest possible expressions of the 
frequency matrix and of the GMFA-GF matrix in the (q, ^-representation enables 
reliable numerical investigation of the consequences coming from the adjustable 
parameters of the model (the degree of hole/electron doping, the energy gap A, the 
hopping parameters). 

Another open question of the GF approach to the solution of the present model is 
the use of the Hubbard operator algebra to get rigorous derivation and simplification of 
the Dyson equation of the complete Green function. As shown previously in [12], the 
self-energy corrections induce a spin fluctuation <i-wave pairing originating in kinematic 
interaction in the second order. 

These investigations are underway and results will be reported in a forthcoming 
paper. 
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